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Q^' We argue that multi-particle production in high energy hadron and nuclear col- 



hsions can be considered as proceeding through the production of gluons in the 
background classical field. In this approach we derive the gluon spectrum immedi- 
^ \ ately after the collision and find that at high energies it is parametrically enhanced 

' by ln(l/x) with respect to the quasi-classical result {x is the Bjorken variable). We 

show that the produced gluon spectrum becomes thermal (in three dimensions) with 
an effective temperature determined by the saturation momentum Qg, T = c Qs/2tt 
during the time ~ 1/T; we estimate c = \/^/2 ~ 1.2. Although this result by 
itself does not imply that the gluon spectrum will remain thermal at later times, it 
has an interesting applications to heavy ion collisions. In particular, we discuss the 
possibility of Bose-Einstein condensation of the produced gluon pairs and estimate 
the viscosity of the produced gluon system. 
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I. INTRODUCTION 

Recently, it was suggested that a fast thermalization in heavy-ion colhsions can occur 
through the gluon radiation off rapidly decelerating nuclei Ij. In that paper two of us 
have pointed out that a pulse of strong chromo-electric field produces Schwinger-like [2I 
radiation with a thermal spectrum. We also discussed an analogy between the Schwinger 
mechanism and the Hawking-Unruh radiation and its application to heavy-ion collisions (see 
also . The macroscopic approach of [H] led to a number of intriguing but qualitative 

results. In the present paper we would like to reconcile the macroscopic approach of (l| with 
the microscopic one based on the QCD parton model. 

The main goal of this paper is to give a picture of the thermalization stage of the process 
of multiparticle production in heavy ion collisions at high energy in the framework of the 
color glass condensate (CGC) approach to high density QCD [g, [Zl, l8|] . The CGC approach 
is based on two principle ideas. The first one is the structure of the parton cascade at high 
energy which is shown in Fig. [H The main contribution to the high energy scattering is 
given by a parton fluctuation in which all partons are strongly ordered in time. Let z be 
the beam direction in the rest frame of the target. The typical lifetime of this fluctuation at 
high energy of the projectile e is large and is proportional to e/m?, where m is the virtuality 
of the fluctuation. In terms of the light-cone variables fc^ = £j ± k^i the life-time of the z-th 
parton is of the order of = x+j = = kf /kf^, where ki± is the transverse momentum 

of the i-th parton. Introducing the rapidity Ui of the parton we can rewrite the lifetime as 

The interaction with the target of the size R destroys the coherence of the parton wave 
function of the projectile. The typical time, which is needed for this, is of the order of 
At ~ i? and is much smaller than the lifetime of all faster partons in the fluctuation: 
At ^ tj. Therefore, this interaction cannot change the momentum distribution of the fast 
parton in the projectile wave function. The influence of the target mostly manifests itself 
in the loss of coherence for majority of the partons; changes in momenta occur only for a 
few very slow ('wee') partons. The 'wee' parton part of the wave function together with 
the interaction with the target could be factorized out while the energy dependence and 
distributions of the fast partons should not depend on the properties of the target. (In the 
following discussion we will assume that the interaction happens at time t = 0.) Such picture 
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FIG. 1: The time structure of the parton cascade for a fast particle (nucleus) in the target rest 
frame. Note, that ordering of the high energy gluons in rapidity y is equivalent to their ordering 
in the space-time rapidity t] = (1/2) ln(x+/x_). In the target rest frame, the target's rapidity 
is r/ = 0, while that of the fast projectile is r] = Y > 0. The relation between the light-cone 
coordinates of the fastest gluon in the cascade is therefore x+ = x_e^^, while for the slowest one 
x+ > x_. 

follows from the parton model and is based on rather general properties of field theories (see 
e. g. ref. jsl); it has been proven in QCD for the BFKL emission 10]. 

The CGC approach adds a very essential new idea to the parton cascade picture. Since 
all partons with rapidity larger than y (see Fig. [1]) live longer than the parton with rapidity 
y, for a dense system such as a nucleus they can be considered as the source of the classical 
field that emits a gluon with rapidity y [oj]. We are going to explore this idea to evaluate the 
parton wave function at the time t = (see Fig. [1]), or to say better just after the interaction, 
when the coherence of the wave function has been destroyed (see section [IT])- Moreover, we 
will argue in Section [TTl that the dominant source of parton production is the longitudinal 
background field; we will also elucidate the origin of this field. 

We then use the same background field approximation to follow the parton system at 
later times. In Section 3 we will argue that the produced parton spectrum assumes the 
three-dimensional thermal form (in a co-moving frame, of course) over the time ~ 1/Qs, 
where Qs is the saturation scale which is a new dimensional parameter that characterizes 
the partonic wave function at t = jo, Q, Q]. We confirm the result of ^ that the effective 
temperature is approximately T = Qs/{2tt). At later times, the partons will interact with 
each other and these interactions finally could create a thermalized system of partons in 
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the true "thermodynamical" sense (for example, with temperature related to the density by 
equation of state), but the consideration of this late kinetic equilibration stage is beyond 
the scope of this paper. We would like to note only that a three-dimensional thermal shape 
of parton distributions should make a true kinetic equilibration easier. 

In this paper we will use also two other key properties of a dense partonic system in 
QCD: 

The first one is the appearance of a new scale (saturation momentum Qg) which 
characterizes the mean transverse momentum of partons in the parton cascade. This mo- 
mentum is proportional to the density of partons (gluons) in the projectile at fixed rapid- 
ity, namely, oc xG{x,Qs)/ttR'^ where xG is the number of gluons with fixed Bjorken 
X = exp(— and R is the transverse size of the projectile. This scale increases with rapidity 
since xG oc (1/a;)^ in the region of low x. It means that the smaller is the value of Qs the 
faster is the parton. Therefore, the parton with rapidity y in Fig. [1] has a mean transverse 
momentum which is much larger than the transverse momentum of all partons moving faster 
than it; thus it can be considered as a probe for the system of fast partons, similar to the 
deep inelastic probe. This observation allows us to consider the production of a parton as 
a process of emission by the frozen system of faster partons; averaging over the quantum 
numbers of incoming hadrons can be done after calculating the cross section. The typical 
configuration of the emitter is such that the transverse sizes are much larger (transverse mo- 
menta are much smaller) than the typical transverse sizes for the emitted parton (transverse 
momentum of emitted parton). 

The second main idea behind the CGC approach is that the quantum emission in each 
stage of the process should give the same result as the emission by the classical field. This 
idea is the cornerstone of the Wilson renormalization group approach in JIMWLK formalism 
(see Ref. |ll|). 



This paper is organized as follows. In Sec. [TTl we explain the origin of the longitudinal 
fields in high energy hadron and nuclear collisions. We then consider the motion of a gluon 
in the external longitudinal color field. The imaginary part of the gluon propagator in an 
external field is related to the cross section for the inclusive gluon production, Eq. lfTTl) . We 
calculate the gluon propagator for an arbitrary external field in Sec. Ill 41 using the WKB 
approximation. In Sec. IIIII we use the derived formulae to calculate the imaginary part of 
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the gluon propagator. Depending on the value of the adiabaticity parameter 7, see (137|) . we 
obtain the gluon spectrum at early times fl50|) and at later times fl5^ . These are the main 
results of our paper. Eq. ( l50l) coincides with the McLerran-Venugopalan formula for gluon 
emission by dense randomly distributed two-dimensional color charges. The corresponding 
saturation scale is given by (H9l) . Eq. (15^ implies that at later times gluon distribution 
is thermal with the temperature determined by the saturation scale (I6OI) . Assuming the 
validity of k± factorization, in Sec. |V] we generalize our formalism to the case of heavy 
ion collisions. In Sec. IVII we consider multiple gluon pair production. Since the gluon 
spectrum at later times is thermal we apply well-known formalism of statistical physics to 
calculate the thermal properties of the produced gluon system. In particular, we observe 
the phenomenon of Bose-Einstein condensation which may solve the long-standing puzzle of 
multiple soft gluon production. We discuss and summarize our results in Sec. I VI II 



II. HIGH ENERGY PARTICLE PRODUCTION BY EXTERNAL FIELDS 



1. Transverse and longitudinal fields of the CGC 

The potential of a charge moving with constant velocity v along the z-axis is given by a 



particular case of the Lienard-Weichert potential (see e. g. |l2l |) 



A, = ^ , '''' (1) 

A± = 0, (2) 

where we introduced the light-cone potential = y4_, A±) with A± = {Aq ± Az)/\/2. 
If the particle is fast, then v ^ 1 and the potential takes form 

r , ^^'"^ (3) 

A_ = = 0. (4) 

The corresponding fields read 

F = ± (1 - ^')^^ 

^ Atx {2x1 + (1 - t;2)xi)3/2 ' ^ ' 

H = V X E. (6) 
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Dirac equation in the background field ©-(jl]) was solved in Ref. [13| with an assumption 
that the fast particle moves freely from a;+ = — oo to = oo. In this case the potentials ([3])- 
(jlj) generate purely transverse mutually orthogonal electric and magnetic fields. The action 
for such a plane wave background vanishes. This implies that there is no pair production in 
a single monochromatic plane wave background j^. 

The initial conditions in our case are different. As explained in the caption of Fig. [H for 
any gluon in the cascade it holds that a;+ > x_ or, equivalently, z > 0. In other words, in the 
target rest frame, all gluons move in the same positive z direction. Therefore, the potential 
exists only in the positive half-plane z > 0. In other words we have to solve the pair 
production problem with the initial condition which explicitly depends on both lightcone 
coordinates x+ and x_. In Sec. Ill 31 we show that such an initial condition generates the 
longitudinal chromoelectric field in addition to the transverse fields mentioned above (see 



( !34|) - (!36l) and below). The existence of longitudinal fie 



has been pointed out previously in Ref. In Ref 



ds in the Color Glass Condensate 
the pair production 



mechanism in heavy-ion collisions by non-perturbative fields has been discussed. 

The longitudinal field Ez is not only generated in a high energy collision, but it gives a 
leading contribution to the pair production amplitude as we are now going to demonstrate. 
Consider a system of fast charges located at coordinates Xj randomly distributed in the 
transverse area of typical size R±. Let us calculate a field created by all these charges at 
the point x^. Assuming for simplicity a continuous distribution of the charge, in the leading 
order in the coupling we have 

f d^xp(x',t^\x-x'\) r (Px'dx' p{^'-.x'l. ^+ + x'- - 2{t~-t^] ) 
^ ^ J 4n \x-x'\ J in x„ - x' + i^iZ^ 

— 2{x-—x'_) 

Typical partons having rapidities y and y' such that y' > y have x_'s satisfying x_ ^ x'_. 
Also, the typical transverse size of a parton decreases down the cascade as x± ~ 1/Qs{y) 
since Qs{y) is an exponentially increasing function of y. Therefore the transverse sizes satisfy 
x'j_ ^ xj_ which implies that the field A{x) does not depend on the transverse size of the 
parton x±: 

^ , , f d^x'^dx'_ v^p(x'_, xl, x+ + x'_ - (|^) 

A Ax) ^ / ^-TTa . (8) 

Eq. ([8]) implies that at high energies the transverse fields experienced by the partons are 
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small compared to the longitudinal ones: 

\E^\ = \H^\ < E,. (9) 

We need to consider the result of (Q with some caution since E± is still enhanced at very 
small values of x_, see (|3]). However in the Lagrangian C = [E"^ — H^)/2 the transverse fields 
indeed give a very small contribution proportional to Ej_ — Hj^ ^ E"^ Ql{y') / Ql{y) <^ E^. 

The pair production probability is proportional to the imaginary part of the effective 
Lagrangian evaluated by considering the quantum fluctuations in the background of the 
external color fields. Therefore, we expect that the pair production will be dominated by 
the longitudinal color fields; we will check this by an explicit calculation below. 



2. Particle production in the background field 

Inclusive production of a gluon with rapidity ?/ in a gluon cascade shown in Fig. [2] can 
be considered as a production of a gluon in a constant background field. Indeed, for this 
gluon all other gluons with larger rapidities are effectively frozen and constitute a constant 
classical field Ez. Therefore, the splitting of a fast gluon into two gluons at rapidity y at 
the time ty {t'y in the complex conjugated amplitude) can be considered as a process of a 
gluon pair production by the field E. As shown in Fig. [2] both gluons propagate in the 
classical background field. Interaction with the target takes much shorter time than the 
gluon emission ty — t'y ^ t — t' . Therefore, the only dynamical role of the interaction with 
the target is to break the coherence of the nuclear wave function and to allow an inclusive 
measurement. This is the reason why we can present the inclusive cross section in a factorized 
form, namely, da/dy oc ipp(ty — ty) -ipTit — t'), where ipp is the probability to find a gluon in 
the projectile. Calculating this probability we could neglect the fact that one gluon interacts 
with the target because of the short interaction time. The unintegrated distribution (fT is 
thus the probability for a gluon to interact with the target. Clearly, this simple factorization 
formula is just another representation for the well-known A;t-factorization formula which 



holds in high density QCD, at least for the interaction o 



a nucleus with a virtual photon or 



19 
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21 
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23|, 



hadron targets, [17] and which has the form (see Refs. 

da AtiNc 1 /" ,,2 /-I A i2\ II 7 \2 \ i-,,^\ 
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In our approach it is convenient to write this formula in a different way 

e^ = j<fkj:\G^2G)lmD{Y-y,p^-k^)lmD{Y-y,k^)My.P±-ki_) (11) 
where 

AttN 

V\G^2G) = a.^^ (12) 

c 

and ImZ)(?/, kt) is the imaginary part of the gluon propagator in the strong classical field. 

Let us consider a target with the transverse size R much smaller than 1/Qs{y) where Qs 
is the saturation momentum. For example, consider the virtual photon target with virtuality 
^ Ql{y)- III this case we can neglect the dependence on kt in ipx in ffTTj) and write 

da 



dy(Pp±_ 

with 



^p{y - y, P±) ^riy, P±) (13) 



ipp{Y-y,p^) = j dH^r\G^2G)lmD{Y-y,p^-k^)lmD{Y-y,p^) (14) 

The dependence on k± is absorbed in the dependence of the classical fields on the trans- 
verse coordinate. In the first approximation we consider the classical fields being independent 
of the transverse coordinate. It means that the gluon propagator is proportional to S^'^\k±) 
and f|T^ can be rewritten as 

MY - y,P±) = T\G ^ 2G) lmD{Y - y,p^) lmD{Y - y,p^) (15) 

The factor l/p^ is included in our definition of D{Y — y,p±)] it must be reproduced for the 
values of p± at which the perturbation theory is valid. 



3. Equation of motion in the background field 



Now we can concentrate our efforts on the calculation of y^p which describes the produc- 
tion of the gluon pair in the strong and constant field. This problem has been investigated 
in detail both in QED and QCD (see review 25| and references therein) and can be solved 
by using the background field method. Let us assume that gluon fields have the structure 



G^ = + W^, 



(16) 
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FIG. 2: Inclusive gluon production at rapidity y in the target rest frame. Crosses mark gluons 
which are on mass shell. ipT and ipp denote the gluon densities for the target and the projectile, 
respectively. The pair of gluons produced at rapidity y at time ty in the amplitude and t'y in the 
complex conjugated one is shown by a bold line. 

where is a classical background field and is a quantum fluctuation. The QCD 
Lagrangian can be expanded around this classical field and it has the following general 



structure [26 1 



Since the second term is equal to zero due to equation of motion for the classical field, our 
Lagrangian has a quadratic form as far as the quantum field dependence is concerned. In 
the case of SU{2) an explicit calculation (see Appendix IVIll) leads to the equation of motion 
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for the quantum field Wn: 



(18) 



where we used the gauge condition D^A^^ = 0. The field configuration discussed in Sec. Ill II 
satisfies this condition since d^A^ = d^A^{x_) = and = 0. 
The tensor for the longitudinal electrical fields is 

/ E,\ 


V -E, / 

where = —d^A^. The components of (ITS!) look as follows 



/i = 0; -{dx - igAxf + igEMz = , 



/i = 3; -(^A - igAxf W^, - igE^W, 



0, 



Introducing W± = Wq ± iW^, we can rewrite (120!) and (12T|) in the form 

-{{dx-igAxf±2gE,)W± = 0. 



(19) 



(20) 
(21) 



(22) 



4. Calculation of a gluon propagator in the background field 



We now turn to solving the Eq. (l22l) . Although A^ is a function of only x_ ( !22|) cannot 
be solved by separation of variables since the initial condition depends on both x+ and x_ 
as has been discussed in Sec. Ill II We can only separate the x± dependence. Thus, we are 
looking for the solution in the form Wn_ = e"*'^"*^-^'^^, where a = ±1. Working in the WKB 



approximation \d-^-Sd-S\ ^ \d-^d-S\ [271, 



we reduce fl22l) to 



2d+Sid-S - gA+{x-)) +pl + 2gaE, = 0, x+ > x. 



(23) 



where d+ = and d- = Eq. (l22l) is a Hamilton- Jacobi equation for motion of 

a charged particle in the background field = A+(x_). The only difference from the 
classical mechanics is the appearance of the spin-dependent term in the right hand side of 
(I23D. 
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In the Hamilton- Jacobi formalism the action S is considered along the true trajectories 
(satisfying Hamilton equations). It is a function of the coordinate x of the final point 
of the trajectory. The action along the true trajectories can be found using the method 
of characteristics. This method was suggested for this class of problems in [3, Q] (for a 



mathematical review see e. g. 30l|). 



A. Solution with a = 

In the spinless case o" = 0, characteristics of Eq. ( l23l) are given by the solution of the 
following set of ordinary differential equations valid at x+ > X- 

dx- 



-2p- , (24) 



dt 

^ = -2{p^ - gA^ix.)) , (25) 

^ = -2p+p_ - 2p_(p+ - gA+) , (26) 
dp- 

dp+ 



0, (27) 

-2p_(7^V(x_), (28) 



dt 

where t is a parameter along the characteristics and we introduced the canonical momenta 
p± as 

P- = -d-S, p+ = -d+S. (29) 

Instead of one of the equations flMl) - fl271) we can use the following equation stemming from 
(123]) and ([29]) 

-2p^{p+ + eA+{x-))+p\ = 0. (30) 

We will use §^ in place of (|27D. 

We can use x_ as a parameter along the characteristics and rewrite (!25l) . (l26l) and (!28|) 
in the following way 



dx+ ^ p+ + gA+{x^) 
dx- p- ' 

^ = 2p^ + gA4x.), (32) 

|± = -,A;(x_). (33) 
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Using Eq. (l30!l the system (l3T]) - (l33|) can be easily integrated with the following result 



p+ = -gA+{x^) + gA+{x+) + pi, (34) 
^_ = ^ / _^^_i^±__^ , (35) 



{pl + gA4x+)y 

S = — qA+{xJ)dx_ + / dx+—7:, , — r . (36) 

7 J ^p% + gA+{x+) ^ ' 

Eq. fl35|l coincides with the equation of motion of a classical test particle of mass p± in the 
external field yl+(x+). In other words, the test particles effectively move under the action 
of the longitudinal electric field = ~A'_^{x+). 

Eq. (1361) gives the action of the test particle along the trajectory (1351) . Its imaginary 
part arises from the pole in the integrand of the second term in the right-hand-side of ( l36l) . 
Integration around the pole in the plain of complex x^ yields the imaginary part. It can 
be calculated replacing the denominator in the first integral in (136!) by Im(p^ + gAj^^)^^ = 
±{in/2)6{pl + gA^) according to the Landau-Cutkosky cutting rule. Additional factor of 
1/2 arises due to the condition a;+ > X-. Define 

T = x^u;, A^ir) = -^f{T), 7 = ^- (37) 

uj gEo 

where is a typical frequency of the external field and E{t = 0) = Eq. With this definitions 
we obtain 

dr 



""^ = '"V 51^7^ = ^Tv^ ■ '''' 

The imaginary part of the action fl38l) corresponds to the pair production. In Fig. [3] we show 
a geometrical interpretation of pair production in the constant background field. 

The physical meaning of the adiabaticity parameter 7 introduced in fl371) is clear: 7 = 
for the static field, while 7 ^ 1 for rapidly oscillating one. Since gEo ~ kf_^ and cu = ki- 
we have the following estimate 

7-f^. (39) 

This estimate for 7 is the quintessence of a quahtative discussions in Sec. [Tll Namely, it 
means that for t = the emission of the gluons is determined by small values of 7 or, in 
other words, by constant electric fields, in which A^{x^) = Eqx+. 
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FIG. 3: Motion of a particle in the constant background field in the light-cone coordinates. Using 

2 

([35|) one can derive the trajectory of the particle = —2^^ (for simplicity we set = 0). 

At X-i- < the particle moves freely along the light cone = until the point 2;_ = X4- = at 
which it tunnels along the line x_ = (Euclidean path, shown by dashed line) to a real trajectory 
at x+ > 0. At this point particles move along x^ > x_ branch of the parabola, while antiparticles 
along the x+ < x_ branch. 

B. a = ±1 case 

In the case of o" = ±1 Eg. (12^ cannot be integrated in general. Equations (ISTIl and ( 132]) 
remain valid in this case. In place of (!32|) we obtain 

^ = -gA'4x_) - ^(p+ - gA4x_))E',{x,) , (40) 

while in place of (!30l) we have 

- 2p4p+ + eA4x^)) + pI + 2g<jE, = . (41) 

Eq. (HUi) can be integrated to yield p^ = ]9+(x_,x+). However, substitution of p+ into 
f lST]) gives an ordinary differential equation which cannot be integrated for an arbitrary field 

y4+(x_). 

We can still investigate the pair production in the two most important cases of constant 
and rapidly decreasing fields. 
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1. If Ez = const then (HOj) reduces to (133|) . Eq. (HTl) then imphes that the solution to (123|) 
with 0" = ±1 is given by fl38l) with shifted transverse momentum Pj_ ^ Pj_ + 2gaEz. 
This can be seen of course directly in (J23l) . 

2. For large a;_ i?^ decreases at least as l/x^J"^ and can be dropped in fj^Oj) and in (jH]) 
bringing us back to the spinless case (l38l) . 

In both cases sum over spins yields an additional factor of 2 in front of fl38l) . 

III. TIME EVOLUTION OF THE CGC WAVE FUNCTION 

1. Model for A+{x-) 

The y4+(x_) potential in ([8]) can in principle be calculated by integrating over the trans- 
verse positions x±i of the highly energetic partons. However, in the present paper we will 
restrict ourselves to a simple model which describes both the small and large x_ behavior 
of the background field. As was discussed in the introduction, at t = = Eq = const 
which implies that A_ = —Eqx^. 

There are two important effects determining the late-time behavior of the chromoelectric 
field. First, the produced gluons start to interact which results in the increase of the gluons' 
k- momentum and hence the field frequency u. We discuss this effect in detail in the 
following subsection. Second, the produced color pairs screen the original color field. The 
invariant mass of the pair provides the mass gap in the excitation spectrum. Therefore, we 
expect the exponential fall-off of the field amplitude. This can be incorporated in a simple 
model 

= :^(l-e— ). (42) 
to 

For this model we can derive the following expression for the imaginary part of the action, 
see (EH]): 

y i^s = J]"^^ , . (43) 

The model potential in fl42p as well as the simple answer of fH51) is, of course, a simpli- 
fication of the real situation. However, it is easy to see that this model incorporates the 
main properties of the parton cascade that we have discussed. In ([8]) the density of the color 
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charge can be approximated by 

p x'^, x+ + x'_ - = c j (fk^e""^^-^'^ 5 {x'_ - uj-^) 5 [kl - Ql) , (44) 

where c is a constant. Note that before the interaction x_ was neghgible since X- ~ 1/^+- 
However, right after the interaction its typical value becomes of the order of x_ ~ uj^^ ~ Q^^ 
which follows from the uncertainty principle Afc_|_ x_ ~ 1 and fl52l) (or fl56l) ). In writing (jH]) 
we also took into account the fact that most of the gluons have transverse momenta of the 
order of ~ Qs- The density given by ( l44l) generates A+(x_) according to ([H]) in the form 

A+(x_)=c j — 2a;^_ + x'l = c x. Ko{x^Qs) , (45) 

where d is another constant. One can see that (H5l) reproduces the main property of the 
model potential of (142|) . Namely , A^{xS) ~ x_ (up to a logarithm) as x_ ^0 and 
A^{x^) ~ exp {—QsX^) as x_ ^ 1. Therefore, we believe that the model potential of P2|) 
reflects the main properties of the structure of the parton cascade in high density QCD 
(CGC). It is worthwhile mentioning that the mass gap turns out to be of the order of the 
saturation momentum and this looks very natural in the CGC approach. 



2. Gluon spectrum at t = 

To calculate the gluon spectrum we have to calculate the imaginary part of the action S 
as explained in Sec. Ill 2[ First, we will calculate the spectrum of produced particles at initial 
time a:+ = and then, in the next section, we will consider later-time particle production. 

It follows from ( l43l) that in the constant field (7 = 0) 

This equation solves the problem of finding the propagator of a gluon with transverse 
momentum px in the strong constant classical field. In the WKB approach we can guarantee 
only the exponential suppression for \m.D{Y — y,p±) and 

lmD{Y-y,p^) oc e"^^-^'"^^ = — e"^ (47) 

as 

Note that the dependence on the spin a canceled out. Substituting (H7|) in ([T^ we obtain 

AttN 

v{Y-y,p^) = Spj;^^e _ (48) 



16 



The coefficient in front of the exponent in (HHIl was chosen based on the physical meaning 
of function Lpp (see Ref. js^). Sp in fHHl) is the transverse area of the projectile and is 
the running QCD coupling. 

Eq. ( HHi) allows us to introduce the saturation scale which is defined to be the mean 
momentum of the produced gluons: 

« = §^- m 

Using this new variable the unintegrated gluon distribution function becomes 

^p{pA.) oc Spj^^ e ^. (50) 



This equation gives the CGC parton density which coincides with the formula suggested 

35I for more detailed calculation of 



by McLerran and Venugopalan in Ref. (see refs. 
CGC parton density), and which illustrates the main property of the CGC approach: the 
entire dependence on rapidity and impact parameter enters only through the saturation 
scale Ql{y,h). 

Therefore, our simple picture leads to the CGC initial condition at t=0. In the next 
section we wish to discuss how the system can develop after losing coherence due to the 
interaction in the final state. 



3. Thermalization by a pulse of the chromoelectric field. 

After losing coherence at t = the fast gluons start to interact 4]- A fast ith gluon in 
the cascade changes its longitudinal momentum and energy according to Newton law 

The second equation states that the energy of a gluon propagating in the constant back- 
ground field, which exists at t = 0, does not change. Eqs. (!5T!) imply that during the time 
of the order of l/Qs the longitudinal momentum changes its value by ~ Qs- This results in 
a variation of both kf and k~ by the same value 

Aki ~ ^kr ~ Q, (52) 

Since kf ^ Qs-, for k~^ it is a small relative change, and can be neglected. However, the 
initially (at t = 0) small value of k~ = k\/kf <^ Qg increases in a strong field up to the 
hard scale Qg. 
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The classical fields will depend on time with the typical frequency of k~ . Therefore, the 
interaction among the fast partons leads to oscillation of the classical fields with a typical 
frequency ~ Q^. However, since the values of kf for the fast partons are still larger than 
Qs we observe that all slow partons (with rapidity y in Fig. [1] and in Fig. [2]) stem from the 
classical emission of the fast partons. 

We now turn to the derivation of the gluon spectrum at later times. It was suggested in 
P] that at later times particles are produced by a pulse of the longitudinal chromoelectric 
field. Indeed, the third equation in (1371) implies that the adiabaticity parameter increases 
with CO. Thus, at later times 7^1. It follows from ( l43l) that in the case of exponentially 
decreasing field (and only in that case) the final spectrum is thermal. The imaginary part 
of the action reads 

where we used = p\/ (2p° ) which is true for the real particles. 
For the imaginary part of the gluon propagator we thus obtain 

lmD{Y-y,p^) = ^e-'^^'^'^ = ^ e"^ , (53) 

The unintegrated gluon distribution is 

^{Y -y,p^) = Spj^^e ^ =^p_^eT. (54) 

Eq. fl5^ implies that at later times the gluon spectrum is thermal with the temperature 



4. Thermalization time 



The initial state of the heavy ion is characterized by the distribution of gluons fl50|) with 
the typical transverse momentum Qg proportional to the strength of the chromoelectric field, 
fH9|) . Since we have assumed that Qs is the only relevant scale, the effective temperature T 
and thermalization time ttherm (over which the spectrum acquires the thermal shape) must 
be related to Qs- To estimate them we will use the following two observations: (i) Due to 
momentum conservation the p+ momentum gained by the emitted particle is equal to the 

momentum lost by the field; (ii) The dominant contribution to the integral of (l38l) comes 
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from times r ~ 7. The value of the adiabaticity parameter 7 ~ 1 marks the transition 
between the gaussian and the thermal distributions. In other words, 

Ap+ = -p% = -UJ, (56) 

r = ujx- ~ 1 (57) 
Eq. (l33ll implies the following estimate 

Ap^ - , (58) 

UJ 

where in the last equation we used fl57|) . Then from fl56l) and (l58l) we estimate the typical 
frequency of the field 

iu = (59) 

Let us now substitute the definitions of the saturation scale (l49ll and the temperature (155|) 
into ( l59l) . The result is 

T^-^Q,. (60) 



The characteristic time over which the field changes is 



and the thermalization time is 



11 

t ~ _ ~ ^ — (61) 



^therm — rp • (^2) 



In the kinematical range of RHIC for the collisions of heavy nuclei Qs — 1 ^ 1-5 GeV |33| 
This translates to T ^ 200 ^ 300 MeV and ttherm =^ 0.6 ^ 1 fm. 



IV. NUCLEAR GLUON DISTRIBUTIONS 



To understand better the approximation that we suggest in this paper we consider here 
the process of deep inelastic scattering off the nuclear target assuming that the nucleus is so 
heavy that we can treat it as a source of the classical field Let us assume that the probe 
is not a virtual photon but is rather a graviton or other particle that can interact with a 
gluon. For such a probe we have two different way of interaction with the target. In the 
first one the probe decays into two gluons and one of them belongs to the classical field of 
the target (see the upper figure in Fig. H]). The second process goes in two steps: the first is 
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Probe 




time 



FIG. 4: The time structure of the deep inelastic scattering in color glass condensate: the upper figure 
describes the quasi-classical McLerran-Venugopalan approach while the lower one corresponds to 
the approach developed in this paper. 

the decay of the probe into two quantum gluons and in the second stage these two gluons 
interact with the classical field as it is shown in low picture in Fig. HI 

For the first process we have McLerran-Venugopalan formula |6| , namely, the distribution 
of produced gluons in the coordinate space looks as 



d?x^dy OLs 



-\xl Ql{x) \n(xl Ql{x)) 



e 4 



(63) 



The second process leads to a formula that at first glance has a quite different form namely 
3ll 

oc ln(l/a;) (l - g-MQ^W M-iQ^W)) (64) 
a-^x±dy \ / 

Eq. (1641) is written in the so-called leading logarithmic log(l/x) approximation (LLA) of per- 
turbative QCD in which we consider only contributions that are proportional to (a^ ln(l/x))" 
such that as ln(l/x) ~ 1 while <^ 1. Since in LLA ln(l/x) ~ l/a^ one may conclude 
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ihat (IMIl gives a contribution of the same order as (l63i) . However, it has been shown in Ref. 



32] that in the saturation region where lxlQl{x) > 1 Eq. ([MD can be rewritten as follows 

oc - In (xl Qlix)) ( 1 - e^i^i "^(^i (65) 

One can see that in the saturation region the contribution of (!65l) which corresponds to 
our approach is parametrically larger than the quasi-classical McLerran-Venugopalan result 
of (l63l) . It is well known that in the wide region of the kinematic variables the mean field 
approximation to Color Glass Condensate leads to the geometrical scaling behavior; namely, 
all experimental obsevables turn out to be functions of one variable C = In (^^l Qs(3^)); in 
which we have the non-linear equation. Even without discussing the exact form of this 
equation one can see that fl63|) is the initial condition for such an equation while fl65l) gives 



its first iteration. The equation itself 



111 ] is based on the idea that each emitted gluon 



with large longitudinal momentum could be treated simultaneously as a quantum and as a 
classical field. Eq. (l65l) is a good illustration of this principle since the quantum emission of 
gluons leads to a result with dN / d?x^dy oc l/ag- 



V. ION-ION COLLISIONS 

For ion- ion collisions we intend to use the kt factorization approach expressed by ( |TOl) . 
This equation has not been proven for CGC; nevertheless, we still think that it provides a 
reasonable starting point, for the following reasons. First, the factorization has been proven 



for large values of transverse momenta 38|] (see also reviews in Ref. 39|). Second, (|TOl) is 
the correct formula for the inclusive production in the case of the BFKL emission (see Ref. 

n 

[19[ and references therein). This fact is very important in understanding why this relation 
could be valid even in the CGC region. Indeed, the BFKL equation has its own, intrinsic 
scale of hardness: the mean transverse momentum of gluons which increases as a function 
of energy. This fact is common for the BFKL and CGC emissions, especially if we recall 
that the BFKL approach is the low parton density limit of the CGC. However, the rigorous 
proof of ffTOj) is still lacking. The theoretical situation as well as physical arguments for such 



factorization have been outlined in Ref. 40|] and we cannot add more at the moment. 
For the ion-ion collision we thus use the following equation 

da da AnNr 



e 



d^p dyd^pt — 1 pi 



T-2 I dk\as'^A{y -y,k]_)LpB{y,{p-k)']_), (66) 
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where is given by (15^ and subscripts A and B refer to the mass numbers of the nuclei. 
This factorization formula can be rewritten in the form of ffTTl) . namely, 



- j d^k^ ImDAiY -y,p^- k^) ImD^(y - y,p^) 

X ImDsiy, p± - k^_) ImDBiy, p±) ■ (67) 



In the first approximation we can integrate over k±_ assuming that the classical fields do not 
depend on the transverse coordinate. Therefore, we have 

da Att A^f, 

lmDA{Y - y,p±) lmDA{Y - y,p±) Im DB{y,p±) lmDB{y,p±) ■ (68) 



dyd'^pi, - 1 

In Fig. [5] we show that the gluon is moving in the fields Ea and Eb in the time interval 
t — t'. In fact, by writing Eq. fl67|) we assumed that the resulting field is just the sum of these 



two fields. It is correct for QED, but not for QCD [6|, [UJ, |35|]. In other words, we assumed 
that during the time interval t — t' both gluons interact with two fields in such a way that 
the resulting propagator is equal to 

lmD{t - t') = ImDA^t - t') luiDBit - f) (69) 

For t = (EHD leads to 

^i^^SASBj2N^^^_ 

dyd^p^ as N^-l ' ^ ' 

the effective saturation scale in ion-ion collisions thus can be inferred from (170|) as 1/Ql = 
^/Q'i,A + expected. 

For t > l/T Eq. §7^ looks differently: 

From (I7T1) we see that we have the same expression as in (jSlj) but with a different tempera- 
ture. Therefore the spectrum is given by 

da _ SaSb tt^N, 
dy(Ppl ~ as 2{N^-r 

with 



p 

" c 



e -^cff , (72) 



111 
— = — + —. (73) 

J eflf J- A J-B 
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FIG. 5: Inclusive production in the collision of the nucleus A with the nucleus B. Crosses mark 
the on-mass- shell gluons. 

VI. STATISTICAL INTERPRETATION OF MULTIPLE PAIR PRODUCTION 

1. Probability of multiple pair production 



It was argued in Ref. 27|, |4l|] that the pair production mechanism allows a statisti- 
cal interpretation. Consider the relative probability of single pair production wi{a,p) = 
exp(— 2ImS'). Assuming that the pairs are produced independently, the absolute probabil- 
ity to produce one pair is then given by 

oo 
n=l 

similar expressions hold for the absolute probabilities to produce n pairs, Wn = Wi{l — 
Yl'^=i'^i)- Let Wo(o", = 1 — '"^r('^'P) be the probability that no pair with quantum 

numbers a, p is produced. The probability conservation condition then reads 

oo oo / _t 

The total probability that the vacuum of a given theory remains unchanged in a given 
volume V during time At is 

Wo = I exp{iC VAt)]"^ = exp(-2 ImC VAt) . (76) 
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On the other hand 27|, |4l|] , 



(77) 



where we used f l75l) . Therefore, 



V 



a,p 



d^p\Yi{l-Wi{a,p)), (78) 



where g = (2(t + l)(A^c — 1) is the degeneracy factor for pairs of particles. The expression 
on the left hand side of flTHj) is nothing but the total production probability in the WKB 
approximation 

1 - Wo ~ 2 ImC VAt . (79) 

When wi is given by thermal distribution ( 153|) the right hand side Eq. ( 1781) is related to the 
thermodynamic potential fip^irs of the produced pairs: 

V 



^pairs B 



(2vr)= 



d pln(l — wi{a,p)) . 



(80) 



Since we work in the approximation in which the background field does not depend on the 
transverse coordinates the particles produced in a given pair are correlated exactly back-to- 
back. Therefore, the thermodynamic potential for single particles Q is just twice the one for 
the pairs. 

2. Bose-Einstein condensation 



For the reasons which will become clear shortly, let us introduce a new notation W/T 
2 ImS. 2n takes the following values at early and later times 



2U/T 



\ V-/T.. 



t ^therm 
t ^therm 



The number of the produced pairs is equal to [42 1 



N 



At=0 



d VT 



g J d^pln{l - e 



{^^-^2IJ)/T^ 



At=0 



(27r)2 J (MIT _ 1 

where we absorbed the additional degeneracy factor 2 in the definition of g. 



(81) 



(82) 
(83) 
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At t = it follows from (IHTj) that the integral in (183|) logarithmically diverges in the 
infrared region in agreement with perturbative QCD. However, the total emitted energy is 
finite 

e = JedN, = ^Ap^QlgEo^am, (84) 

where we have restricted ourselves to the central rapidity region where e = p± = v^P- = 
\/2p+, pz = 0. We can estimate Ap^ — gEot. Also during the longitudinal expansion ~ t. 
Therefore, during the early stages after the collision the energy flows from the field to the 
soft particles as S ~ t^. 

The gluon number becomes finite as soon as t > 0. Indeed, when t changes from to 
^therm it follows from flHTl) that n = d\nW/dliap± decreases from 2 to 1. We have 

where A varies from at t = to T at t = ttherm and 

The distribution in ( !86l) at n = 1 has the form of a Bose-Einstein distribution with a 
vanishing chemical potential, yU = 0. We thus expect the Bose-Einstein condensation of 
gluons to occur at temperatures lower than the critical temperature Tq. For the sake of 
simplicity let us assume that t is close to ttherm so that ti > 1 and A > T. To calculate Tq 
let us note that (l85ll cannot be used for counting the number of particles which carry zero 



42| . where To is defined as 



transverse momentum = at T < Tq 

3(27r)^Ar \ 
VQgEtTi 

The number of particles with zero momentum (in the condensate) equals 

2 



N{p^ = {]) = N 



(88) 



lierm- 



whence is the total (finite) number of particles. 

The critical temperature decreases with time. Let us estimate its value at t = tthe 
Using (H9l) and ( l62l) we get (y'-Eottherm ~ 1- Assume V = Sxt where Sa is the transverse cross 
sectional area of the nucleus (for simplicity we assume a central collision). Then 

12N\ 1/2 



To < To (ttherm) = ( ) " (89) 
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The total number of hadrons produced at ?/ = at RHIC is about ~ 1000. Using 
Sa = 7r(7 fm)^ and g = 2 ■ 3 • 8 we obtain To(ttherm) ~ 260 MeV. Therefore, just after the 
system is thermahzed, a significant fraction of gluons may form a Bose-Einstein condensate. 

The Bose-Einstein condensation of soft gluons in high energy QCD leads to a remarkable 
consequence. Recall that the typical correlation length inside the high energy hadron is 
rather small Ac — 1/Qs- This implies that the gluon emission with long wavelengths A = 
> -^c is suppressed because it decouples from the hadron wave function, similarly to the 
decoupling of a large wavelength signal from a small antenna. Therefore, one is led to predict 
a deficit of soft gluons at high energies, in a stark contradiction with the experimental data. 
The phenomenon of Bose-Einstein condensation solves this puzzle since it allows piling up 
of soft gluons. 



3. Viscosity of the parton system 

We have argued in Sec. IIII 31 and Sec. IIII4I that at t > ttherm the produced partons have 
3D thermal distributions with an effective temperature T. The number of produced particles 
per unit volume is large n ~ l/a^ since they were part of the classical fields in the initial 
wave functions before the collision. This observation is an important argument in support 



of the hydrodynamical description of the parton system at later times [43 1. 

The typical transverse momentum of a parton is {p±) ~ T. Recall that the temperature 
T is proportional to the saturation scale Qs{y) which is an exponential function of rapidity. 
Therefore, the temperature varies with rapidity. As a consequence, the average value of 
transverse momentum {p±) significantly varies between different rapidity layers. The dif- 
ference in the transverse momentum distributions along the longitudinal axis of rapidity 
amounts to the shear viscosity^. 

The shear viscosity can be estimated as (we keep only parametric dependence while 
omitting all numerical factors) 

^ = {P±)\-—, (90) 
n na 

where a ~ otg/Ql is the scattering cross section for a parton in the classical background 



^ We would like to thank Ben Svetitsky for bringing our attention to this consequence of our approach. 
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field. The number of particles per unit volume is 

where ~ 1/Qs is the longitudinal extent of the system. Using Ql ~ agxG/SA we then 
estimate 

- ~ 1. (92) 

n 

This estimate implies the parametric smallness of viscosity which comes about as a conse- 
quence of high occupation number of gluons in the initial wave function. In contrast, in 



pQCD the shear viscosity is parametrically enhanced rj/n l/a1 44 1. 



VII. DISCUSSION AND SUMMARY OF THE RESULTS 

In this paper we have developed an approach to particle production based on the principle 
idea of CGC: the gluon with a rapidity can be considered as emitted by the classical fields 
that are composed of all faster partons with y > i/q. We showed that in such an approach 
the gluons at the moment of collision are emitted by the classical longitudinal fields {E^), 
which are created by fast particles during very short time ~ 1/e after the collision {e is the 
particle energy in the laboratory frame). We found the relation (H9l) between the momentum 
scale of dense partonic system Qg and the strength of the classical field E^. The inclusive 
distribution at t = is given by fl50l) and turns out to be the same as has been expected in 
the CGC approach (McLerran-Venugopalan model) j^. 

At later times t > 1/ ^/2tt Qg, one has to consider the time dependence of the classical 
fields. We followed through the evolution of the system assuming that the main source 
of the produced gluons is still the classical field created by faster partons. It turns out 
that the momentum distribution of the produced gluons has a three-dimensional thermal 
spectrum given by (15^ with T = (l/4:\/27r)Qs{y) for the collision of two identical nuclei 
at midrapidity. Therefore, the CGC approach led to a thermal spectrum of emitted gluons 
with an effective temperature which depends on the rapidity of emitted gluons. 



It was argued in Ref. 45| that the perturbative dynamics may not be adequate for the 
description of the late-time processes in a high-energy heavy-ion collisions as it does not 
lead to the thermalization as anticipated on general grounds. In the present paper we 
circumvent that result by suggesting a non-perturbative mechanism of thermalization. The 



27 



non-perturbative nature of the obtained results can be clearly seen in Eq. (l38i) which exhibits 
non-analytic dependence on the coupling g. 



The dependence o 



(see for example 
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temperature on rapidity may trigger instability of the gluon system 
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481] and references therein) and speed up thermalization process. 



Perhaps at late times the instability driven thermalization can compete with pair production 
by strong fields discussed in this paper. This problem warrants further investigation. 

Another problem left beyond the scope of the present paper is understanding at what 
time the hydrodynamic description becomes valid. It seems reasonable to asssume that for 
times later than ttherm ^ 1/T we could apply the viscous hydrodynamic description. Indeed, 
we showed that for these times we have a 3D thermal distributions in each slice of rapidity 
which is a pre-condition for using the hydrodynamic approach. On the other hand, the 
average transverse momenta {pt) — T are quite different in the two neighboring slices in 
rapidity due to the dependence of T on rapidity. Therefore, we can expect a considerable 
difference in parton momentum distributions in different rapidity slices which amounts to 
viscosity. We have argued that the CGC initial conditions lead to the parametrically small 
shear viscosity rj ~ 0{1) as opposed to the perturbative result, rj ~ 0{l/a1). It should be 
mentioned that matching the CGC energy-momentum tensor with that of an almost perfect 
fluid yielded similar results |49|]. 

Although we performed our calculations for the SU{2) gauge theory, we believe that all 
the qualitative features of the derived results will remain valid for the realistic color group 
SU (3) as well. Calculations of the pair production effect in a constant chromoelectric field 
of SU{3) have been recently done in Ref. 50|, ImI]- Unlike SU{2) there are two Casimir 
operators in SU (3) which yield a more complicated dependence of the pair production effect 
on E. 

A new related general approach to particle production in field theories coupled to strong 
external sources has been recently formulated in Ref. 52| where the particular example 
of X(f)^ theory has been discussed. It may yield new insights into the problem of particle 
production problem in QCD as well. 

It is interesting to note that calculation of inclusive e'^e' production in QED can be done 
in exactly the same way as was followed to calculate the gluon production in this paper. 
Indeed, a fast moving system in QED is characterized by large transverse fields which lead to 
bremsstrahlung production of photons which is a classical process. There is also production 
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of e^e~ pairs which is a typical quantum process. The QED variant of the CGC approach 
states that at high energies the inclusive production is dominated by the emission of e^e" 
pairs in the classical photon field and not by the quantum emission of virtual photons. 
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APPENDIX A: EQUATION OF MOTION OF A VECTOR PARTICLE IN AN 

EXTERNAL FIELD IN SU{2) 

Let = be the background classical field. We are looking for the equations of motion 
of the vector particle PV^ = {A^^ + iA^^/\/2 in the background field A^. The Lagrangian is 

1=1 

Using identity 

1=1,2 

[id.Al - d^Al) + g{AlAl - AlADf-r [{d,Al - 9.^) - g{AlAl - A^^)]' , (A2) 
and expanding the 3-component of the strength tensor 

{Fl^f = (d^K - d.Alf + 2g{d^Al - d.Al) 2t W; + 0{Wi) (A3) 

we obtain 

£ = -\id^A, - d,A^)^ - ig{d^A, - d,A^)w;W, -\\D^W, - D^W^W (A4) 
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where = + igA^. The corresponding equation of motion is 

[-D^V + D^D, + 2tgF^,]W'' = 0. (A5) 
Assuming D^W^ = we have 

dlW, - 2igA^d^W, - ig{d^A^)W, - g^ AlW, + 2gEaW, = 0, (A6) 
which is equivalent to ( fTSjl . 
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